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Abstract 

We report an ab initio-hased study of mesoscopic quantum transport in chemically doped 
graphene nanoribbons with a width up to 10 nm. The occurrence of quasibound states related to 
boron impurities results in mobility gaps as large as 1 eV, driven by strong electron-hole asym- 
metrical backscattering phenomena. This phenomenon opens new ways to overcome current 
limitations of graphene-based devices through the fabrication of chemically-doped graphene 
nanoribbons with sizes within the reach of conventional lithography. 

The basic principle of CMOS field effect transistor (CMOS-FET) is to drive a high current 
density in its ON-state, whereas electrostatic gating blocks almost completely the current flow in 
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the so-called OFF-state. The current ratio between OFF and ON states is one of the key parameters 
determining the device performance.- This electrostatic control of the CMOS-FET is strongly 
efficient owing to the use of semiconducting materials, which present a sufficiently large electronic 
bandgap to electrostatically deplete or accumulate charges in the conducting channel. Recently, 
the synthesis of low dimensional carbon based materials has opened alternatives to silicon-based 
electronics devices. Indeed, by graphite exfoliation^ or by sublimation of SiC,-^ the fabrication of 
single layer two-dimensional (2D) graphene or stacked graphene multilayers was achieved, with 
exceptionally high values of reported charge mobilities (> 100.000 cm^V^^s^^) close to the Dirac 
point.-^^ Notwithstanding, with 2D graphene being a zero-gap semiconductor,^ a poor field effect 
efficiency with small ON/OFF current ratio was obtained.-^*^ 

E-beam lithography and oxygen plasma etching allow to fabricate graphene nanoribbons (GNRs),-^»^ 
down to a few tens of nanometers in width, opening the possibility for bandgap engineering through 
electronic confinement. However, current limitations of state-of-the-art lithographic techniques 
hardly allow to fabricate GNRs with sufficiently large energy bandgaps.— ^ ^^ ^ ^ To further re- 
duce the GNRs lateral size down to a few nanometers, alternative chemical approaches have been 
developed and competitive device performances have been reported.—*^ These approaches remain 
however hardly compatible with massive integration of interconnected devices at the wafer scale. 

Chemical doping allows to develop new applications such as chemical sensors or more ef- 
ficient electrochemical switches.—*^ The conductivity modulation due to chemical doping and 
functionalization of 2D graphene has been explored numerically ,-22a^ without showing significant 
improvement of device performances. Several theoretical studies have also investigated the effect 
of doping on the properties of single-layer GNRs. The impact of a single boron impurity on the 
electronic properties of GNRs was analyzed by first-principles methods in Refs.— ^^^i Other works 
have considered the GNRs bandgap tuning by doping or functionalization,'^^^'^ but only gaps of 
a few tens of meV have been predicted. Theoretically, two types of GNRs with highly symmetric 
edges have been described, namely, zigzag (zGNRs) and armchair (aGNRs). To date, most efforts 
have focused on the spin-transport properties of zGNRs, while the possibilities for charge transport 
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of doped aGNRs have somehow been neglected. Besides, due to the computational cost of ab initio 
calculations, the width of the GNRs studied in those works is in general well below 3 nm. 

In this Article, we report on a quantum transport study of chemically doped GNRs up to 10 nm 
in width. By combining first-principles calculations with tight-binding models, the effect of chemi- 
cal doping on charge conduction is explored for ribbons with a length up to the micron and random 
distribution of substitutional boron impurities. The joint effects of backscattering together with the 
enhanced localization due to low dimensionality are shown to produce strong electron-hole trans- 
port asymmetry and large mobility gaps, even in the low doping limit. These chemically modified 
graphene-based materials, within reach of current lithographic techniques, should therefore allow 
for a strongly enhanced current density modulation under electrostatic gating. 

We focus on the case of aGNRs. This type of ribbon symmetry has theoretically been proposed 
as energetically more stable and more resistant to edge defects when compared to zGNRs,— i^"^ 
which seems to be confirmed by experimental observations.— As reported by ab initio calcula- 
tions,— >^^^ aGNRs are found to be always semiconducting with an energy gap that varies up 
to 0.5 eV for ribbons of similar width, depending on the exact number of dimer chains within 
the unit cell. Unlike zGNRs, the ground state of aGNRs is not spin-polarized, p-type (n-type) 
doping can be achieved by the incorporation of boron (nitrogen) atoms in substitution within the 
sp^ -type carbon matrix, as reported in Ref.-^^ for carbon nanotubes and very recently in Ref.^^ in 
the case of 2D graphene. First-principles calculations of the effect of a single boron impurity on 
the electronic properties of aGNRs have been reported in Ref.— for aGNRs with widths between 
2.3 and 4.2 nm. Following Ref.,^^ we refer to a aGNR with N dimers contained in its unit cell 
as an N-aGNR. Both types of possible electronic structures, namely the pseudo-metallic ribbons 
20-aGNR and 35-aGNR (i.e., ribbons predicted to be metallic within a simple nearest-neighbor 
tight-binding model), and the semiconducting 34-aGNR were addressed. The modification of the 
electronic structure of the aGNRs by a single boron (or nitrogen) impurity turns out to be strongly 
dependent on the position of the dopant with respect to the ribbon edge, displaying a significant 
electron-hole asymmetry. -^^'^^ This phenomenon results from the interplay between wavefunction 
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symmetries and screening effects that are system-dependent. Regardless of their exact position 
across the ribbon width, boron impurities are found to preserve their well-known acceptor char- 
acter (or hole-doping ability) in carbon, with vanishingly weak impact of the transport properties 
in the first conduction band close to the Fermi level (E^), which almost preserves conductance 
quantization. 

The large variation of resonant energies with the acceptor (donor) chemical impurity position 
allows however to infer that a random distribution of impurities will likely result in an uniform 
reduction of conductance over the occupied (unoccupied) states part of the first conduction plateau, 
yielding marked electron-hole asymmetry and charge mobility gaps intrinsic to p-type (n-type) 
doping. The present work is the first theoretical study based on ab initio atomistic models providing 
confirmation to such prediction. 

Due to the large size of our unit cells, a statistical study of the mesoscopic transport using the 
fully ab initio Hamiltonian~ is not possible here. We have used instead a combination of first- 
principles calculations with tight-binding models, that has been satisfactorily employed before for 
doped carbon nanotubes.^^ As a first step, the scattering potential induced by the dopant has been 
extracted from those first-principles density functional theory (DFT) calculations. The analysis of 
the ah initio onsite and hopping terms around the dopant, as obtained from self-consistent DFT 
calculations within an atomic-like basis, allows to build a simple nearest-neighbor tight-binding 
(nn-TB) model. The substitutional impurity modifies mainly the onsite terms, creating a potential 
well on a typical length scale of about 10 Angstrom (see Ref.— and Ref.^*^ for a detailed analysis). 
A modification of the C-C TB hopping at the edges is also included in the TB Hamiltonian follow- 
ing Ref.— to simulate the edge passivation, that yields a small energy gap at the charge neutrality 



point (CNP). As shown in [Figure T] the conductance for a single dopant (in various ribbon symme- 
tries and impurity location) computed by our TB models nicely reproduces the ab initio transport 
calculations performed with the TABLIER code^*^ for an extensive set of dopant positions and 
ribbons width. The conductance of long GNRs is then calculated within the Landauer-Biittiker 
formalism. Order(N) decimation techniques are used to explore the various transport properties 
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Figure 1 : Conductance profile as a function of energy for the single-dopant case, both within the 
full ab initio approach (dashed red line) and the tight-binding model (solid black line) for the (a) 
20-, (b) 35-, (c) 31- and (d) 26-aGNRs; the dopant has been placed at the (a) center, (b) midway 
between center and edge, (c) first neighbor to atom at edge, and (d) edge positions, respectively. 

for randomly doped ribbons as long as 1 /im. By combining the accuracy of first-principles to de- 
scribe the local electronic structure with the numerical efficiency of a nn-TB model, our transport 
methodology thus allows for a realistic exploration of the mesoscopic transport in randomly doped 
GNRs. 

We first analyze the impact of different doping rates on ribbons of about 4 nm in width, namely 
the pseudo-metallic 35-aGNR and the semiconducting 34-aGNR. Ribbons with a length up to 1 /im 
are considered, and impurities are uniformly distributed over the whole ribbon length and width. 



with a restriction preventing the overlap of the scattering potentials of individual dopants. Figure 2 
shows the conductance as a function of energy for the (a) 35- and (b) 34-aGNRs, for doping rates 
between ^ 0.02% and 0.2%. In this formalism, electrodes are treated as semi-infinite, perfect 
GNRs, and disorder is included only in the region between the electrodes (channel). Since the 
scattering potential induces an overall shift of the CNP within the channel, the CNP of electrodes 
and channel are not aligned, leading to a displacement of the doping-induced gap with respect to 
that of the perfect ribbon. To properly analyze our results we have thus simply proceeded to a rigid 
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Figure 2: (a) Average conductance as a function of energy for the pseudo-metallic 35-aGNR for 
doping rates ^ 0.02%, 0.05%, 0.1% and 0.2% (from top to bottom). The dashed black line cor- 
responds to the ideal (undoped) case. The averages have been performed over ^ 500 disorder 
realizations with a ribbon length of ^ 1 jim. The inset shows the unit cell of the 35-aGNR with 
passivating H atoms. N indicates the number of dimer chains in a N-aGNR. (b) Same as in (a, 
main frame) for the semiconducting 34-aGNR. Inset: Schematic view of a randomly doped GNR. 

shift of the conductance profile of the doped system to superimpose it to that of the ideal case. 
As a result of the acceptor-like character of the impurity states induced by the boron dopant, the 
conductance is affected in a clear asymmetric fashion for energy values below or above the CNP,— 
as evidenced by the opening of a large mobility gap that extends well beyond the first conductance 
plateau below the small initial electronic bandgap. The mobility gap width reaches almost 1 eV, 
of the order of the silicon energy gap. 

Interestingly, despite the smaller size of the initial gap of the pseudo-metallic 35-aGNR with 
respect to that of the 34-aGNR, the magnitude of the doping-induced mobility gap in both rib- 
bons is very similar for the same ribbon length and doping rate. This fact, as well as the higher 
conductance values for electron transport in the 35-aGNR (due to the robustness against scattering 
processes of the linear antibonding n* band in pseudo-metallic ribbons), make them even more 
suitable for their use in electronic devices. 

Since the width of these GNRs is reaching the capability limits of standard lithographic tech- 
niques, we further explore the case of aGNRs with a width of ^ 10 nm, namely the 80- and the 
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Figure 3: Main panel: Same as Figure 2 for the semiconducting 81-aGNR and three selected 
doping rates 0.02%, 0.05% and 0.2%, from top to bottom). Inset: Same as in main frame for 
the pseudo-metallic 80-aGNR. 

81-aGNR. For these systems, a doping rate below 0.2% is enough to achieve a mobility gap of ^ 



1 eV (see Figure 3 while keeping high conductance values in the conduction band for both types 



of wide GNRs. In this case, higher doping rates are required to achieve the full suppression of 
conductance in the valence band. This is due to the fact that backscattering in aGNRs is highly 
dependent on the dopant position, with a maximum impact when the dopant is close to (or right at) 
the ribbon edge.— As a consequence of a uniform random distribution of dopants, the probability 
to find an impurity close to the edges decreases with increasing ribbon width, and will be, for a 
80-aGNR, almost one half than for a 35-aGNR, thus leading to a lower impact on the conductance 
for a similar doping rate. 

The asymmetry in the electron vs. hole conduction— is also spectacularly evidenced by 



scrutinizing at the energy-dependent transport regimes. Figure 4 (main frame) shows the length 
dependence of the conductance for selected individual randomly doped 35-aGNRs, with a doping 
rate of ~ 0.05%, at a energy = 0.25 eV below (a) and above (b) the CNR For energy values lying 
in the conduction band (b), the conductance for various random configurations is seen to slowly 
decay with ribbon length with values staying close to the conductance quantum (Go = 2e^/h), 
pinpointing a quasi-ballistic regime. Though all 1 -dimensional systems show localization in the 
presence of disorder, for lower doping rates quantum interference effects are too weak to allow 
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Figure 4: Main frame: Conductance as a function of length for a set of doping realizations of the 
35-aGNR at 0.25 eV below (a) and above (b) the CNP for a doping rate of ^ 0.05%. Curves with 
same color in (a) and (b) correspond to the same position of the dopants along the ribbon. Inset: 
Estimated localization lengths as a function of doping rate for the 34- and 35-GNRs (at E = -0.2 
eV), and the 80- and 81-aGNRs (at E = -0.1 eV). Solid lines are fittings to the calculated values. 



localization to manifest itself in the length scale considered in this study. In contrast, for energies 
lying in the valence band, a strong exponential decrease of the conductance is found (a), with large 
fluctuations associated to different defects positions. The exponential decay of the conductance 
with the length of the ribbon is related to the Anderson localization, that has already been observed 
experimentally at room temperature in defected metallic carbon nanotubes,— due to their long 
phase coherence lengths. 

In the inset of Figure 4[ the estimated localization lengths {^) for hole transport at different 
doping rates are reported. A statistical average over more than 500 disorder samples is performed 
to extract the ^ from (InG/ Go) = —^/L, with L the ribbon length. The localization length ranges 
within 10-300 nm, depending on the ribbon width and doping density, and is inversely proportional 
to the doping density. ^ is also found to further increase with the ribbon width but in a non-linear 
fashion. This comes from the non- uniformity of the disorder potential with the dopant position. 

In the case of carbon nanotubes, the localization length is predicted to increase linearly with 
the size of the system,— a trend confirmed experimentally.^-^ A similar scaling law was derived in 
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pseudo-metallic GNRs.— This can be explained by invoking a generalization of Thouless relation 
for quasi-one dimensional systems, ~ N±le, where A^^ is the number of conducting modes avail- 
able at a certain energy and 4 is the elastic mean free path. In the first plateau around the CNP, 
ribbons of different widths present the same number of channels {N± = 1), and Ig is predicted to 
increase with the number of atoms in the unit cell,— thus leading to an increase of ^ . However, 
this is true only in the presence of a uniform disorder potential, which is not the case here. We note 
that similar power law scaling of <^ has been observed for the case of edge disorder.— 

In conclusion, we have performed a numerical transport study of boron doped armchair GNRs 
with widths up to ^ 10 nm, and lengths reaching the micron scale and doping rates from 0.02% to 
0.2%. Depending on the energy of charge carriers, the transport can vary from quasiballistic to a 
strongly localized regime, as a result of a strong electron/hole asymmetry induced by the chemical 
doping. This phenomenon leads to mobility gaps in the order of 1 eV, while the conductance in 
the conduction bands remains large, with values close to those of the ideal ribbon for low enough 
doping rates. This opens an unprecedented way to improve the performances of graphene-based 
devices since the obtained asymmetry should also manifest in a considerably higher ON/OFF cur- 
rent ratio. One notes that electron-hole conduction asymmetry in graphene devices by chemical 
functionalization has been recently reported by Farmer and coworkers,— whereas Raman studies 
of graphene-FETs allow for efficient monitoring of doping effects.— 
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